Evidence of natural selection and dominance of SARS-CoV-2 variant Lambda (C.37) over variants of concern in Cusco, Peru.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage C.37 (Lambda) has spread rapidly in Peru and other Latin American countries. However, most studies in Peru have focused on Lima, the capital city, without knowing the dynamics of the spread of the variant in other departments. Cusco, Peru, is one of the most popular departments in the country for tourists, so the introduction of new variants of SARS-CoV-2 might occur despite closure of the borders. Therefore, in this work, we analyzed the variants circulating in Cusco. The aim of this work was to better understand the distribution of SARS-CoV-2 lineages circulating in Cusco and to characterize the genomes of these strains. To this end, 46 SARS-CoV-2 genomes from vaccinated and unvaccinated patients were sequenced in the first half of 2021. The genomes were analyzed using phylogenetic and natural selection methods. Phylogenetic trees from Cusco showed dominance of the Lambda lineage over the variants of concern (VOCs), and there was no clustering of variants by district. Natural selection analysis revealed mutations, mainly in the spike protein, at positions 75, 246, 247, 707, 769, and 1020. In addition, we found that unvaccinated patients accumulated more new mutations than did vaccinated patients, and these included the F101Y mutation in ORF7a, E419A in NSP3, a deletion in S (21,618-22,501), and a deletion in ORF3a (25,437-26,122).

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations.

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Pandemic-scale phylogenomics reveals the SARS-CoV-2 recombination landscape.

Accurate and timely detection of recombinant lineages is crucial for interpreting genetic variation, reconstructing epidemic spread, identifying selection and variants of interest, and accurately performing phylogenetic analyses1-4. During the SARS-CoV-2 pandemic, genomic data generation has exceeded the capacities of existing analysis platforms, thereby crippling real-time analysis of viral evolution5. Here, we use a new phylogenomic method to search a nearly comprehensive SARS-CoV-2 phylogeny for recombinant lineages. In a 1.6 million sample tree from May 2021, we identify 589 recombination events, which indicate that around 2.7% of sequenced SARS-CoV-2 genomes have detectable recombinant ancestry. Recombination breakpoints are inferred to occur disproportionately in the 3' portion of the genome that contains the spike protein. Our results highlight the need for timely analyses of recombination for pinpointing the emergence of recombinant lineages with the potential to increase transmissibility or virulence of the virus. We anticipate that this approach will empower comprehensive real-time tracking of viral recombination during the SARS-CoV-2 pandemic and beyond.

Distinct evolutionary trajectories of SARS-CoV-2-interacting proteins in bats and primates identify important host determinants of COVID-19.

The coronavirus disease 19 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a coronavirus that spilled over from the bat reservoir. Despite numerous clinical trials and vaccines, the burden remains immense, and the host determinants of SARS-CoV-2 susceptibility and COVID-19 severity remain largely unknown. Signatures of positive selection detected by comparative functional genetic analyses in primate and bat genomes can uncover important and specific adaptations that occurred at virus-host interfaces. We performed high-throughput evolutionary analyses of 334 SARS-CoV-2-interacting proteins to identify SARS-CoV adaptive loci and uncover functional differences between modern humans, primates, and bats. Using DGINN (Detection of Genetic INNovation), we identified 38 bat and 81 primate proteins with marks of positive selection. Seventeen genes, including the ACE2 receptor, present adaptive marks in both mammalian orders, suggesting common virus-host interfaces and past epidemics of coronaviruses shaping their genomes. Yet, 84 genes presented distinct adaptations in bats and primates. Notably, residues involved in ubiquitination and phosphorylation of the inflammatory RIPK1 have rapidly evolved in bats but not primates, suggesting different inflammation regulation versus humans. Furthermore, we discovered residues with typical virus-host arms race marks in primates, such as in the entry factor TMPRSS2 or the autophagy adaptor FYCO1, pointing to host-specific in vivo interfaces that may be drug targets. Finally, we found that FYCO1 sites under adaptation in primates are those associated with severe COVID-19, supporting their importance in pathogenesis and replication. Overall, we identified adaptations involved in SARS-CoV-2 infection in bats and primates, enlightening modern genetic determinants of virus susceptibility and severity.

Evidence for selection on SARS-CoV-2 RNA translation revealed by the evolutionary dynamics of mutations in UTRs and CDSs.

RNA translation is the rate-limiting step when cells synthesize proteins. Elevating translation efficiency (TE) is intuitively beneficial. Particularly, when viruses invade host cells, how to compete with endogenous RNAs for efficient translation is a major issue to be resolved. We collected millions of worldwide SARS-CoV-2 sequences during the past year and traced the dynamics of allele frequency of every mutation. We defined adaptive and deleterious mutations according to the rise and fall of their frequencies along time. For 5'UTR and synonymous mutations in SARS-CoV-2, the selection on TE is evident near start codons. Adaptive mutations generally decrease GC content while deleterious mutations increase GC content. This trend fades away with increasing distance to start codons. Mutations decreasing GC content near start codons would unravel the complex RNA structure and facilitate translation initiation, which are beneficial to SARS-CoV-2, and vice versa. During this evolutionary arms race between human and virus, SARS-CoV-2 tries to improve its cis elements to compete with host RNAs for rapid translation.

Computational Inference of Selection Underlying the Evolution of the Novel Coronavirus, Severe Acute Respiratory Syndrome Coronavirus 2.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that recently emerged in China is thought to have a bat origin, as its closest known relative (BatCoV RaTG13) was described previously in horseshoe bats. We analyzed the selective events that accompanied the divergence of SARS-CoV-2 from BatCoV RaTG13. To this end, we applied a population genetics-phylogenetics approach, which leverages within-population variation and divergence from an outgroup. Results indicated that most sites in the viral open reading frames (ORFs) evolved under conditions of strong to moderate purifying selection. The most highly constrained sequences corresponded to some nonstructural proteins (nsps) and to the M protein. Conversely, nsp1 and accessory ORFs, particularly ORF8, had a nonnegligible proportion of codons evolving under conditions of very weak purifying selection or close to selective neutrality. Overall, limited evidence of positive selection was detected. The 6 bona fide positively selected sites were located in the N protein, in ORF8, and in nsp1. A signal of positive selection was also detected in the receptor-binding motif (RBM) of the spike protein but most likely resulted from a recombination event that involved the BatCoV RaTG13 sequence. In line with previous data, we suggest that the common ancestor of SARS-CoV-2 and BatCoV RaTG13 encoded/encodes an RBM similar to that observed in SARS-CoV-2 itself and in some pangolin viruses. It is presently unknown whether the common ancestor still exists and, if so, which animals it infects. Our data, however, indicate that divergence of SARS-CoV-2 from BatCoV RaTG13 was accompanied by limited episodes of positive selection, suggesting that the common ancestor of the two viruses was poised for human infection.IMPORTANCE Coronaviruses are dangerous zoonotic pathogens; in the last 2 decades, three coronaviruses have crossed the species barrier and caused human epidemics. One of these is the recently emerged SARS-CoV-2. We investigated how, since its divergence from a closely related bat virus, natural selection shaped the genome of SARS-CoV-2. We found that distinct coding regions in the SARS-CoV-2 genome evolved under conditions of different degrees of constraint and are consequently more or less prone to tolerate amino acid substitutions. In practical terms, the level of constraint provides indications about which proteins/protein regions are better suited as possible targets for the development of antivirals or vaccines. We also detected limited signals of positive selection in three viral ORFs. However, we warn that, in the absence of knowledge about the chain of events that determined the human spillover, these signals should not be necessarily interpreted as evidence of an adaptation to our species.

Genomic characterization of the dominating Beta, V2 variant carrying vaccinated (Oxford-AstraZeneca) and nonvaccinated COVID-19 patient samples in Bangladesh: A metagenomics and whole-genome approach.

Bangladesh is experiencing a second wave of COVID-19 since March 2021, despite the nationwide vaccination drive with ChAdOx1 (Oxford-AstraZeneca) vaccine from early February 2021. Here, we characterized 19 nasopharyngeal swab (NPS) samples from COVID-19 suspect patients using genomic and metagenomic approaches. Screening for SARS-CoV-2 by reverse transcriptase polymerase chain reaction and metagenomic sequencing revealed 17 samples of COVID-19 positive (vaccinated = 10, nonvaccinated = 7) and 2 samples of COVID-19 negative. We did not find any significant correlation between associated factors including vaccination status, age or sex of the patients, diversity or abundance of the coinfected organisms/pathogens, and the abundance of SARS-CoV-2. Though the first wave of the pandemic was dominated by clade 20B, Beta, V2 (South African variant) dominated the second wave (January 2021 to May 2021), while the third wave (May 2021 to September 2021) was responsible for Delta variants of the epidemic in Bangladesh including both vaccinated and unvaccinated infections. Noteworthily, the receptor binding domain (RBD) region of S protein of all the isolates harbored similar substitutions including K417N, E484K, and N501Y that signify the Beta, while D614G, D215G, D80A, A67V, L18F, and A701V substitutions were commonly found in the non-RBD region of Spike proteins. ORF7b and ORF3a genes underwent a positive selection (dN/dS ratio 1.77 and 1.24, respectively), while the overall S protein of the Bangladeshi SARS-CoV-2 isolates underwent negative selection pressure (dN/dS = 0.621). Furthermore, we found different bacterial coinfections like Streptococcus agalactiae, Neisseria meningitidis, Elizabethkingia anophelis, Stenotrophomonas maltophilia, Klebsiella pneumoniae, and Pseudomonas plecoglossicida, expressing a number of antibiotic resistance genes such as tetA and tetM. Overall, this approach provides valuable insights on the SARS-CoV-2 genomes and microbiome composition from both vaccinated and nonvaccinated patients in Bangladesh.

Influenza viruses and coronaviruses: Knowns, unknowns, and common research challenges.

The development of safe and effective vaccines in a record time after the emergence of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a remarkable achievement, partly based on the experience gained from multiple viral outbreaks in the past decades. However, the Coronavirus Disease 2019 (COVID-19) crisis also revealed weaknesses in the global pandemic response and large gaps that remain in our knowledge of the biology of coronaviruses (CoVs) and influenza viruses, the 2 major respiratory viruses with pandemic potential. Here, we review current knowns and unknowns of influenza viruses and CoVs, and we highlight common research challenges they pose in 3 areas: the mechanisms of viral emergence and adaptation to humans, the physiological and molecular determinants of disease severity, and the development of control strategies. We outline multidisciplinary approaches and technological innovations that need to be harnessed in order to improve preparedeness to the next pandemic.

Enhanced High Mutation Rate and Natural Selection to Produce Attenuated Viral Vaccine with CRISPR Toolkit in RNA Viruses especially SARS-CoV-2.

The best and most effective way to combat pandemics is to use effective vaccines and live attenuated vaccines are among the most effective vaccines. However, one of the major problems is the length of time it takes to get the attenuated vaccines. Today, the CRISPR toolkit (Clustered Regularly Inerspaced Short Palindromic Repeats) has made it possible to make changes with high efficiency and speed. Using this toolkit to make point mutations on the RNA virus's genome in a coculture of permissive and nonpermissive cells and under controlled conditions can accelerate changes in the genome and accelerate natural selection to obtain live attenuated vaccines.

Comparative genomic analysis reveals varying levels of mammalian adaptation to coronavirus infections.

Severe acute respiratory coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is of zoonotic origin. Evolutionary analyses assessing whether coronaviruses similar to SARS-CoV-2 infected ancestral species of modern-day animal hosts could be useful in identifying additional reservoirs of potentially dangerous coronaviruses. We reasoned that if a clade of species has been repeatedly exposed to a virus, then their proteins relevant for viral entry may exhibit adaptations that affect host susceptibility or response. We perform comparative analyses across the mammalian phylogeny of angiotensin-converting enzyme 2 (ACE2), the cellular receptor for SARS-CoV-2, in order to uncover evidence for selection acting at its binding interface with the SARS-CoV-2 spike protein. We uncover that in rodents there is evidence for adaptive amino acid substitutions at positions comprising the ACE2-spike interaction interface, whereas the variation within ACE2 proteins in primates and some other mammalian clades is not consistent with evolutionary adaptations. We also analyze aminopeptidase N (APN), the receptor for the human coronavirus 229E, a virus that causes the common cold, and find evidence for adaptation in primates. Altogether, our results suggest that the rodent and primate lineages may have had ancient exposures to viruses similar to SARS-CoV-2 and HCoV-229E, respectively.

COVID-19 pandemic and "survival of the fittest".

Occurrence and recurrence of COVID-19 cases have been observed globally. The complex relationship of host-pathogen and the environment plays a vital role in understanding the widespread recurrence of the SARS-CoV-2 among humans. Though the pathobiology of the disease is not completely understood, it is well established that COVID-19 poses a greater threat to individuals with co-morbidities and a weakened immune system. The article deals with the notion of innate immunity, natural selection, and the survival of the fittest during the COVID-19 outbreak. The article also attempts to introduce the concept of "lifestyle and cultural immunity" that needs to be addressed and incorporated at an early stage of childhood to boost up the human immune system. The communication further discusses the role of vaccination and micro-organisms pre-existing in the environment which are required to enhance the immunity of an individual.

Nucleocapsid mutations R203K/G204R increase the infectivity, fitness, and virulence of SARS-CoV-2.

Previous work found that the co-occurring mutations R203K/G204R on the SARS-CoV-2 nucleocapsid (N) protein are increasing in frequency among emerging variants of concern or interest. Through a combination of in silico analyses, this study demonstrates that R203K/G204R are adaptive, while large-scale phylogenetic analyses indicate that R203K/G204R associate with the emergence of the high-transmissibility SARS-CoV-2 lineage B.1.1.7. Competition experiments suggest that the 203K/204R variants possess a replication advantage over the preceding R203/G204 variants, possibly related to ribonucleocapsid (RNP) assembly. Moreover, the 203K/204R virus shows increased infectivity in human lung cells and hamsters. Accordingly, we observe a positive association between increased COVID-19 severity and sample frequency of 203K/204R. Our work suggests that the 203K/204R mutations contribute to the increased transmission and virulence of select SARS-CoV-2 variants. In addition to mutations in the spike protein, mutations in the nucleocapsid protein are important for viral spreading during the pandemic.

Circulation and Evolution of SARS-CoV-2 in India: Let the Data Speak.

The COVID-19 pandemic is a global challenge that impacted 200+ countries. India ranks in the second and third positions in terms of number of reported cases and deaths. Being a populous country with densely packed cities, SARS-CoV-2 spread exponentially. India sequenced ≈0.14% isolates from confirmed cases for pandemic surveillance and contributed ≈1.58% of complete genomes sequenced globally. This study was designed to map the circulating lineage diversity and to understand the evolution of SARS-CoV-2 in India using comparative genomics and population genetics approaches. Despite varied sequencing coverage across Indian States and Union Territories, isolates belonging to variants of concern (VoC) and variants of interest (VoI) circulated, persisted, and diversified during the first seventeen months of the pandemic. Delta and Kappa lineages emerged in India and spread globally. The phylogenetic tree shows lineage-wise monophyletic clusters of VoCs/VoIs and diversified tree topologies for non-VoC/VoI lineages designated as 'Others' in this study. Evolutionary dynamics analyses substantiate a lack of spatio-temporal clustering, which is indicative of multiple global and local introductions. Sites under positive selection and significant variations in spike protein corroborate with the constellation of mutations to be monitored for VoC/VoI as well as substitutions that are characteristic of functions with implications in virus-host interactions, differential glycosylation, immune evasion, and escape from neutralization.

Insights on the SARS-CoV-2 genome variability: the lesson learned in Brazil and its impacts on the future of pandemics.

Since the beginning of the SARS-CoV-2 spread in Brazil, few studies have been published analysing the variability of viral genome. Herein, we described the dynamic of SARS-CoV-2 strains circulating in Brazil from May to September 2020, to better understand viral changes that may affect the ongoing pandemic. Our data demonstrate that some of the mutations identified are currently observed in variants of interest and variants of concern, and emphasize the importance of studying previous periods in order to comprehend the emergence of new variants. From 720 SARS-CoV-2 genome sequences, we found few sites under positive selection pressure, such as the D614G (98.5 %) in the spike, that has replaced the old variant; the V1167F in the spike (41 %), identified in the P.2 variant that emerged from Brazil during the period of analysis; and I292T (39 %) in the N protein. There were a few alterations in the UTRs, which was expected, however, our data suggest that the emergence of new variants was not influenced by mutations in UTR regions, since it maintained its conformational structure in most analysed sequences. In phylogenetic analysis, the spread of SARS-CoV-2 from the large urban centres to the countryside during these months could be explained by the flexibilization of social isolation measures and also could be associated with possible new waves of infection. These results allow a better understanding of SARS-CoV-2 strains that have circulated in Brazil, and thus, with relevant infomation, provide the potential viral changes that may have affected and/or contributed to the current and future scenario of the COVID-19 pandemic.

SARS-CoV-2 evolution in animals suggests mechanisms for rapid variant selection.

SARS-CoV-2 spillback from humans into domestic and wild animals has been well documented, and an accumulating number of studies illustrate that human-to-animal transmission is widespread in cats, mink, deer, and other species. Experimental inoculations of cats, mink, and ferrets have perpetuated transmission cycles. We sequenced full genomes of Vero cell-expanded SARS-CoV-2 inoculum and viruses recovered from cats (n = 6), dogs (n = 3), hamsters (n = 3), and a ferret (n = 1) following experimental exposure. Five nonsynonymous changes relative to the USA-WA1/2020 prototype strain were near fixation in the stock used for inoculation but had reverted to wild-type sequences at these sites in dogs, cats, and hamsters within 1- to 3-d postexposure. A total of 14 emergent variants (six in nonstructural genes, six in spike, and one each in orf8 and nucleocapsid) were detected in viruses recovered from animals. This included substitutions in spike residues H69, N501, and D614, which also vary in human lineages of concern. Even though a live virus was not cultured from dogs, substitutions in replicase genes were detected in amplified sequences. The rapid selection of SARS-CoV-2 variants in vitro and in vivo reveals residues with functional significance during host switching. These observations also illustrate the potential for spillback from animal hosts to accelerate the evolution of new viral lineages, findings of particular concern for dogs and cats living in households with COVID-19 patients. More generally, this glimpse into viral host switching reveals the unrealized rapidity and plasticity of viral evolution in experimental animal model systems.

Influence of treatment with neutralizing monoclonal antibodies on the SARS-CoV-2 nasopharyngeal load and quasispecies.

OBJECTIVES: To evaluate the impact of neutralizing monoclonal antibody (mAb) treatment and to determine whether the selective pressure of mAbs could facilitate the proliferation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with spike protein mutations that might attenuate mAb effectiveness. PATIENTS AND METHODS: We evaluated the impact of mAbs on the nasopharyngeal (NP) viral load and virus quasispecies of mAb-treated patients using single-molecule real-time sequencing. The mAbs used were: Bamlanivimab alone (four patients), Bamlanivimab/Etesevimab (23 patients) and Casirivimab/Imdevimab (five patients). RESULTS: The NP SARS-CoV-2 viral load of mAb-treated patients decreased from 8.2 log10 copies/mL before administration to 4.3 log10 copies/mL 7 days after administration. Five immunocompromised patients given Bamlanivimab/Etesevimab were found to have mAb activity-reducing spike mutations. Two patients harboured SARS-CoV-2 variants with a Q493R spike mutation 7 days after administration, as did a third patient 14 days after administration. The fourth patient harboured a variant with a Q493K spike mutation 7 days post-treatment, and the fifth patient had a variant with a E484K spike mutation on day 21. The emergence of the spike mutation was accompanied by stabilization or rebound of the NP viral load in three of five patients. CONCLUSION: Two-mAb therapy can drive the selection of resistant SARS-CoV-2 variants in immunocompromised patients. Patients given mAbs should be closely monitored and measures to limit virus spread should be reinforced.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide.

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Ancient Adaptative Evolution of ACE2 in East Asians.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been posing an unprecedented challenge to global public health. SARS-CoV-2 and several other coronaviruses utilize angiotensin-converting enzyme 2 (ACE2) as their entry receptors. The ACE2 gene has been found to experience episodic positive selection across mammals. However, much remains unknown about how the ACE2 gene evolved in human populations. Here, we use population genetics approaches to investigate the evolution of the ACE2 gene in 26 human populations sampled globally. We find the ACE2 gene exhibits an extremely low nucleotide diversity in the East Asian populations. Strong signals of selective sweep are detected in the East Asian populations, but not in the other human populations. The selective sweep in ACE2 is estimated to begin in East Asian populations ∼23,600 years ago. Our study provides novel insights into the evolution of the ACE2 gene within human populations.